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ABSTRACT 

We present measurements of the pairwise velocity dispersion (PVD) for different 
classes of galaxies in the Sloan Digital Sky Survey (SDSS). For a sample of about 
200,000 galaxies with redshifts in the interval 0.01 < z < 0.3, and r-band magnitudes 
Mo.i r between —16 and —23, we study the dependence of the PVD on galaxy proper- 
ties such as luminosity, stellar mass (M*), colour (g — r), 4000 A break strength (D4000), 
concentration index (C), and stellar surface mass density (/x*). The luminosity depen- 
dence of the PVD is in good agreement with the results of Jing & Borner. for the 
2dFGRS catalog, once the photometric bandpass difference between the two surveys 
is taken into account. The value of a\% measured at k = Ih Mpc -1 decreases as a 
function of increasing galaxy luminosity for galaxies fainter than L* , before increasing 
again for the most luminous galaxies in our sample. Each of the galaxy subsamples 
selected according to luminosity or stellar mass is divided into two further subsam- 
ples according to colour, D4000, C and /i*. We find that galaxies with redder colours 
and higher D4000, C, and /z* values have larger PVDs on all scales and at all lumi- 
nosities/stellar masses. The dependence of the PVD on parameters related to recent 
star formation (colour, D4000) is stronger than on parameters related to galaxy struc- 
ture (C, /x*), especially on small scales and for faint galaxies. The reddest galaxies 
and galaxies with high surface mass densities and intermediate concentrations have 
the highest pairwise peculiar velocities, i.e. these move in the strongest gravitational 
fields. We conclude that the faint red population located in rich clusters is responsible 
for the high PVD values that are measured for low-luminosity galaxies on small scales. 

Key words: galaxies: clusters: general - galaxies: distances and redshifts - cosmology: 
theory - dark matter - large-scale structure of Universe. 



1 INTRODUCTION 

The pairwise peculiar velocity dispersion (PVD) measures 
the relative motions of galaxies, and so reflects the action of 
the local gravitational fields. The PVD is thus an important 
quantity for probing the mean mass density fio of the Uni- 
verse and the clustering power on small scales. It is widely 
used to constrain cosmogonic models and galaxy formation 
recipes (Davis et al. 1985). As a well-defined statistical quan- 
tity, the PVD can be estimated from redshift surveys, either 
by modeling redshift distortions in the two-point correlation 
function (2PCF), or by measuring the redshift-space power 
spectrum (Davis & Peebles 1983; Jing & Borner 2001b). 
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The first method relies on the fact that the peculiar mo- 
tions of galaxies affect only their radial distances in redshift 
space. Thus the information for peculiar velocities along the 
line of sight can be recovered by modeling the redshift-space 
2PCF £ {s) (r p , 7r) as a convolution of the real-space 2PCF 
£(r) with the distribution function of the pairwise velocity 
/(«lfl): 

£ W (r p ,Jr) = J f{v 12 )t(^rl + {-K-vr2) 2 ) dv 12 , (1) 

where v\2 = vwij'p,^) is the pairwise peculiar velocity; r v 
and 7T are the separations perpendicular and parallel to the 
line of sight. The real-space correlation function £(r) is usu- 
ally inferred from the projected 2PCF w p (r p ), which is a 
simple Abel transform of £(r). However, the form of /(U12) 
is not known from a rigorous theory. Based on observational 
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(Davis & Peebles 1983; Fisher et al. 1994) and theoretical 
considerations (e.g. Diaferio & Geller 1996; Sheth 1996), an 
exponential form is usually adopted: 
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where «i2 is the mean and a\ 2 is the dispersion of the one- 
dimensional peculiar velocicites. Assuming an infall model 
for vi2(r), the pairwise velocity dispersion ai2 can then be 
estimated as a function of projected separation r p by com- 
paring the observed ^ 3 \r p ,n) with the modeled one. The 
form of the infall model usually adopted is based on the 
self-similar solution and is a good approximation to the real 
infall pattern in CDM models with [3 = ag^o' 6 ~ 0.5 (Jing, 
Mo & Borner 1998, hereafter JMB98). Using this method 
and various redshift samples of galaxies, the measurement 
of PVD has been carried out by many authors (Davis & Pee- 
bles 1983; Mo, Jing & Borner 1993; Fisher et al. 1994; Zurek 
et al. 1994; Marzke et al. 1995; Somerville, Davis & Primack 
1997). There have been significant variations in the results 
of these studies. In many cases the measurement errors were 
underestimated and there were also problems due to "cosmic 
variance" as a result of the limited volume contained in the 
early redshift surveys (Mo, Jing & Borner 1997). As a result, 
the early determinations are of limited power in constrain- 
ing large-scale stucture theories. The first accurate determi- 
nation of the PVD (ui2 = 570 ± 80kms _1 at a projected 
separation r p = 1 h^ 1 Mpc) was presented by JMB98 using 
the Las Campanas redshift Survey (LCRS), and has been 
verified by Zehavi et al. (2002) with the early data release 
of the Sloan Digital Sky Survey (SDSS) and by Hawkins et 
al. (2003) with the two-degree field galaxy redshift survey 
(2dFGRS). This value is substantially higher than the ear- 
lier results based on smaller surveys (340±40kms^ 1 ) given 
by Davis & Peebles (1983). The large difference in these re- 
sults is due to the fact that the value of the PVD is very 
sensitive to the presence (or absence) of rich clusters in a 
sample (Mo, Jing & Borner 1993), and so a sample as large 
as LCRS is needed in order to give reliable estimates. 

The PVD of galaxies can also be determined from the 
redshift-space power spectrum. Although the power spec- 
trum is just the Fourier transform of the 2PCF, it is some- 
times advantageous to work with the power spectrum in- 
stead of with the 2PCF. The advantage of using the redshift- 
space power spectrum to determine the PVD is that it is sim- 
ple and accurate to model the infall effect (Jing & Borner 
2004, hereafter JB04). It is also not necessary to assume a 
functional form for the real space P(k), as is usually required 
for the 2PCF. The relation between the power spectrum in 
redshift space (k, /i) and that in real space P{k) can be 
written as (Peacock & Dodds 1994; Cole, Fisher & Weinberg 
1995): 

P (s) {k,fi) = P(k)(l + /3^) 2 D(kfia 12 (k)). (3) 

Here k is the wavenumber, /i the cosine of the angle between 
the wavevector and the line of sight, and f3 the linear redshift 
distortion parameter. The first multiplies after the power 
spectrum on the right hand side of Eq.(3) is the Kaiser lin- 
ear compression effect (Kaiser 1987), and the term D is the 
damping effect caused by the random motion of the galaxies. 
The damping function D, which should generally depend on 
k, fi and <7i2(fc), was found by Jing & Borner (2001a) to be 



a function of one variable k/j,ai2(k) only. They also found 
that, although the functional form of D[kfiai2(k)] depends 
on the cosmological model and bias recipes, for small k (large 
scales) where D > 0.1, D[k^,ai2{k)] can be approximately 
expressed by the Lorentz form, 



D(k^iG 12 {k)) = 



k 2 ii 2 ai2(k) 2 



(4) 



We will use equations (3) and (4) to derive a 12 and P(k) 
from a measurement of the redshift-space power spectrum 
P' s '(fc, ji). However, when applying this method, there are 
a couple of important points that should be kept in mind. 
First, the Lorentz form is the Fourier transform of the expo- 
nential form of /(U12) when <ri2 is a constant and v\2 is zero. 
This means that ai2(k) drived in the Fourier space is gen- 
erally not the same as ai2(r) defined in the configuration 
space. Furthermore, even if ai2 is a constant, the redshift 
power spectrum of equation (3) may lead to an unphysical 
distribution of pairwise velocities in the configuration space 
when the infall is not zero (i.e. /3 / 0) (Scoccimarro 2004). 
Nevertheless, based on extensive tests of numerical simula- 
tions, Jing & Borner (2001a) demonstrated that, o\2(r) and 
<Ti2(fc) have similar dependences on scale, and are different 
only by 15 per cent if r = 1/k is used in the comparison. 
Therefore, the quantity (T12 (fc) can still be regarded as a good 
indicator for the pairwise velocity dispersion. Second, it is 
usually difficult to determine o\2 and (3 simultaneously, be- 
cause there exists a strong degeneracy in determining these 
parameters (Peacock et al. 2001) from pW(i, //) at small 
scales. Moreover, there could be some luminosity depen- 
dence in [3. JB04 investigated this by using the luminosity 
dependence of the bias parameter b given in (Norberg et 
al. 2002) , and showed that their determinations of the PVD 
are robust to reasonable changes of the /3 values. Thus a 
viable way to determine the PVD from P^(k,fi) is to fix 
f3 at a reasonable estimate and then determine P(k) and 
o\2 from the data. In this way, the real-space power spec- 
trum P(k) is mainly determined from P ( -"\k,^L = 0), and 
the measured PVD <ri2 is usually a function of k. Adopt- 
ing Eq. (4) for the damping function and setting ji = 0.45, 
Jing & Borner (2001b) measured the PVD <7i2(fc) for the 
LCRS and found that the measured PVD is consistent with 
the result reported by JMB98 from the redshift space 2PCF 
measurement. 

Both of the methods described above for determining 
the PVD of galaxies are based on the measurement of galaxy 
clustering. Previous studies of galaxy clustering in the local 
Universe have established that it depends on a variety of 
quantities, including luminosity (e.g. Norberg et al. 2001; Ze- 
havi et al. 2002, 2005), colour (e.g. Zehavi et al. 2002, 2005), 
morphology (e.g. Zehavi et al. 2002; Goto et al. 2003), and 
spectral type (e.g. Norberg et al. 2002; Budavari et al. 2003; 
Madgwick et al. 2003). It is thus expected that there should 
be some dependence of the galaxy PVD on these physical 
properties. Using the 2dFGRS and the method based on 
the redshift-space power spectrum, JB04 present the first 
determination of the PVD for galaxies in different luminos- 
ity intervals. The analysis leads to a surprising discovery: the 
relative velocities of the faint galaxies are very high, around 
700 km s -1 , reaching values similar to those found for the 
brightest galaxies in the sample. The relative velocities at 
intermediate luminosities M* — 1 (M* is the characteristic 
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luminosity of the Schechter [1976] function) exhibit a well- 
defined minimum near 400 km s -1 . This discovery indicates 
that faint galaxies, and the brightest ones, are both prefer- 
entially in massive haloes of galaxy cluster size, but most 
M* galaxies are in galactic scale haloes. Such a luminosity 
dependence of the PVD is not reproduced by the halo oc- 
cupation model of Yang et al. (2003), although this model 
reproduces well the luminosity dependence of the clustering. 

In this paper, we use the second data release of the 
SDSS to examine the dependence of the PVD not only on lu- 
minosity, but also, for the first time, on other physical prop- 
erties such as the stellar mass, star formation rate, stellar 
surface mass density and concentration. These dependences 
will provide useful clues about how the luminosity depen- 
dence of the PVD is produced in the real Universe. As we 
will see, the high PVD values measured at faint luminosities 
can be explained by a population of red galaxies located in 
rich groups and clusters, which is otherwise sub-dominant in 
many other clustering measures such as the 2PCF and mea- 
surements of galaxy-galaxy lensing. The results are consis- 
tent with the conjecture made by JB04 for this population. 
In a separate paper (Li et al. 2006, hereafter Paper I), we 
have measured the 2PCF for the galaxy samples studied in 
this paper. Taken together, our determinations of the 2PCF 
and the PVD should be extremely useful for constraining 
models of galaxy formation. 

In the following section we briefly describe the mea- 
surements of the 2PCFs of galaxies that are presented in 
Paper I. In §3 we describe the method and the results for 
determining the PVD from the redshift power spectrum. We 
summarize our results in the final section. 

Throughout this paper, we assume a cosmological 
model with the density parameter Qq = 0.3 and the cosmo- 
logical constant Ao = 0.7. To avoid the — 51og 10 h constant, 
the Hubble's costant h = 1, in units of 100 kms _1 Mpc _1 , 
is assumed throughout this paper when computing absolute 
magnitudes. In this paper, the quantities with a superscript 
asterisk are those at the characteristic luminosity /mass (e.g. 
characteristic luminosity L*), while the quantities with a 
subscript asterisk are those of stars (e.g. stellar mass M„). 



2 MEASUREMENTS OF THE TWO-POINT 
CORRELATION FUNCTIONS 

In Paper I, we described in detail our procedure for con- 
structing the observational and random samples and for 
measuring the two-point correlation functions in redshift 
space. The samples analyzed here are exactly the same as 
those presented in Paper I and the reader is referred to this 
paper for more details. 

The redshift-space 2PCF £^(r p ,7r) for each subsample 
is measured using the Hamilton (1993) estimator, 



4:DD(r p ,TY)RR(r p ,TY) 



- 1. 



(5) 



Here DD(r p ,ir) is the count of data-data pairs with per- 
pendicular separations in the bins r p ± 0.5Ar p and with 
radial separations in the bins ir ± 0.5A7T, RR(r p ,ir) and 
DR(r p ,ir) are similar counts of random-random and data- 
random pairs, respectively. When computing 2PCFs, some 
possible biases such as the variance in mass-to-light ratio 



and fibre collisions have been carefully corrected (see Paper 
I for a detailed description). Figs.l and 2 show the contours 
of ^ s \r p , n) for galaxies of different luminosities and of dif- 
ferent stellar masses respectively. Both the effect of redshift 
space distortions on small scales (often called the Finger-of- 
God effect) and the infall effect (Kaiser 1987) on large scales 
are clearly visible: on small scales ^ s \r p ,ir) is stretched in 
the 7r-direction and on large scales the contours are squashed 
along the line-of-sight direction. 

In Paper I, we presented the projected two-point cor- 
relation function w p (r p ) which is estimated from ^ s '(r p , ty). 
With this estimator, we investigated the clustering proper- 
ties of galaxies with various physical properties. In the next 
section, we will present their real-space power spectra and 
pairwise velocity dispersions. 



3 THE POWER SPECTRUM AND THE 
PAIRWISE VELOCITY DISPERSION 

3.1 Method 

Following the method of JB04, we obtain for each subsample 
the redshift-space power spectrum as 

P^(k,n) = 27vJ2^A7v l rl ] A\nr p , j e s) (r p , J ,7v l ) 
cos (k^m ) Jo (k p r p j) W g (r p j , m ) , 

where Jo is the zeroth-order Bessel function (Jing & Borner 
2001b); k p and k n are the wavenumbers perpendicular and 
parallel to the line of sight, and are related to k and /x with 



k — \/ k p H- k-n , — k^jk. 



(7) 



Here 7T; runs from —40 to 40 h~ x Mpc with Am = 1 h' 1 Mpc 
and r p j from 0.1 to 50 ft -1 Mpc with Alnr Pij = 0.23. W g 
is the Gaussian window function, which is used to improve 
the measurement by down- weighting £' s ' (r p , n) at the larger 
scales and has the form 



W g (r p ,Tr) = exp 



2S* 2 



(8) 



The smoothing scale is set to be S = 20 bT 1 Mpc following 
JB04. With these parameters, JB04 showed that the real 
space P(k) and the PVD a"i2(fc) can be reliably measured 
at scales > 0.1 ft Mpc -1 and k > 0.2 h Mpc~* respectively. 

From the measurement of the redshift-space power spec- 
trum P (s) (k, fi) , we can determine the real space power spec- 
trum P(k) and the PVD ai2 simultaneously by modeling 
the measured P^ s \k,iJ,) using Eq.(3). As discussed in §1, it 
is difficult to simultaneously derive /3, ai2 and P(k) from 
p( s \k,n). We therefore fix (3 = 0.45 as a reasonable esti- 
mate (Tegmark et al. 2004) and determine P(k) and a\i 
from the data. We will show in the next section that the de- 
terminations of the PVD are robust to reasonable changes 
of the (3 values. 

3.2 The dependence on luminosity 

In Fig. 3 the symbols with error bars show the measured 
power spectrum in redshift space P^ (k, fx). The four pan- 
els correspond to four different luminosity intervals (sam- 
ples L5, L7, L9, and Lll in Table 1 of Paper I). In each 
panel, the values of k range from 0.1 h Mpc -1 at the top, to 
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-17.0<M„, r <-16.0 -18.0<M , r <-17.0 - 19.0<M„, r <- 18.0 




-20.0<M„, r <-19.0 -21.0<M„, r <-20.0 -22.0<M„, r < -21 .0 
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" 1 " 





-20 20 -20 20 -20 20 

r p (h _1 Mpc) r p (h _1 Mpc) r p (h"'Mpc) 

Figure 1. The rcdshift space two-point correlation function £' s ) (r p ,ir) in different luminosity intervals, as indicated. The contour levels 
are increased by factors of 2 from the lowest (0.1875) to the highest (48.0). 



8.5<log 10 (M„/M o )< 9.0 9.0<log 10 (M + /M o )<9.5 9.5<log 10 (M + /M o )< 10.0 




r p (h~'Mpc) r p (h _1 Mpc) r p (h _1 Mpc) 



Figure 2. The redshift space two-point correlation function £' s ' (r^, tt) in different stellar mass intervals, as indicated. The contour levels 
are increased by factors of 2 from the lowest (0.1875) to the highest (48.0). 
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Figure 3. The redshift space power spectrum p( s )(fc,/z) (symbols with error bars) for galaxies in different luminosity intervals, as 
indicated above each panel. The lines are the best fits of Eq.(3) to the data. In each panel, from top to bottom, the wavenumbcr k is 
0.1, 0.16, 0.25, 0.4, 0.63, 1.0, 1.58, 2.51, and 3.98 h Mpc" 1 , respectively. 



4.0 h Mpc -1 at the bottom, with an increment of A log 10 k — 

0. 2. The error bars plotted in this figure and in all sub- 
sequent figures are estimated by the bootstrap resampling 
technique (Barrow, Bhavsar, & Sonoda 1984). We generate 
100 bootstrap samples and compute the power spectrum 
for each sample using the weighting scheme (but not the ap- 
proximate formula) given by Mo, Jing, & Borner (1992). The 
errors are the scatter of the power spectra among these boot- 
strap samples. The solid lines are the best fits obtained by 
applying Eq.(3) to the data. The real-space power spectrum 
P(k), determined from the modelling with Eq.(3), is dis- 
played in Fig. 4 for 6 luminosity samples (Samples L5, L7, L9, 
L10, Lll and L13 of Table 1 in Paper I). To see the system- 
atic change with the luminosity, a function P(k) — (60/fe) 1 ' 4 
is plotted as the dashed line in each panel. From both figures, 
it is seen that there exists a luminosity dependence of the 
real-space power spectrum: P^ s \k, = 0) and P(k) increase 
with increasing luminosity on all scales, with the luminos- 
ity dependence becoming stronger for bright galaxies and on 
small scales. This result is consistent with that in JB04 for 
the 2dFGRS, and with our analysis of the 2PCF in Paper 

1. Such a dependence on luminosity was also found for the 
2PCF by Norberg et al. (2002) using the 2dFGRS and by 
Zehavi et al. (2005) using the SDSS. The real-space power 
spectrum P(k) is approximately a power law for the range of 
k considered here. The dotted lines in Figure 4 are the best 
power-law fits to the measured P(k) with the best-fitting 
parameters indicated in each panel. P(k) decreases with k 
approximately as k~ 1A , with the slopes of the brightest sam- 
ples being somewhat shallower (~ fc -1,2 ). Moreover, Figure 
4 shows that there is a slight change in the slope of P(k) at 
k w 1 h Mpc -1 for faint galaxies and at k ~ 0.5 h Mpc -1 for 



bright ones. Such a change in slope is much less pronounced 
in JB04, but it is evident here due to the improvement in the 
observational data. As pointed out by JB04, such a change 
can be understood as the transition between the scales where 
the pair counts are dominated by galaxy pairs in the same 
halo to those where galaxy pairs are mostly in separate ha- 
los. This result may indicate that fainter galaxies tend to 
reside in smaller halos, while brighter galaxies are found in 
massive halos, but a more detailed analysis is needed to ex- 
plain the subtle changes of the power spectrum shape with 
luminosity. It is interesting to note that such a transition was 
also found by Zehavi et al. (2004) in the projected two-point 
correlation function of galaxies. 

Fig. 5 shows the PVD o"i2(fc), which is measured simul- 
taneously with P(k), for galaxies in different luminosity in- 
tervals. The errors of 0"i2(fc) are also estimated by bootstrap 
resampling technique. It can be seen that for the fc-values 
used here, <7i2(fc) is a well-determined quantity. As discussed 
in §1, there could be a luminosity dependence of the lin- 
ear redshift distortion parameter /3, so a fixed value of (3 
might have some impact on our results. We investigate this 
by using the luminosity dependence of the bias parameter 
b(M)/b* = 0.895 + 0.150(L/L*) - 0.040(M - AT) given in 
Tegmark et al. (2004) and (3* = 0.45 (quantities with an 
asterisk are defined to be at the characteristic luminosity 
L*). The results are plotted in Fig. 5 and agree with those 
obtained for /? = 0.45 within the error bars. The largest 
differences are seen for bright galaxies on very large scales. 
This indicates that our cri2(fc) measurements are robust to 
reasonable changes of the (3 values. 

The dependence of the PVD on luminosity is clearly 
seen in this figure. The observed PVD for L* galaxies (the 
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Figure 4. The real space power spectrum P(k) for galaxies in different luminosity intervals, as indicated. The dotted lines are the 
power-law fits with the best-fitting parameters indicated in each panel. To guide the eye, the power spectrum P(k) = (60/fc) 1 ' 4 is plotted 
as dashed line in every panel. 
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Figure 5. The pairwise velocity dispersion ai2 for galaxies in different luminosity intervals, as indicated. The symbols with errorbars are 
for results that obtained by fixing the redshift distortion parameter f3 = 0.45, while the lines are for those with /3 varying with luminosity 
as in Tegmark et al. (2004). 
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Figure 6. The pairwise velocity dispersion <ri2 as a function of 
luminosity, measured at k = 0.25,0.5,1.0 and 4. Oh Mpc" 1 . For 
clarity, small constants were added to the Mo.i r values for the 
curves of k = 0.25,0.5 and 4.0/i Mpc -1 . 
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Figure 7. The pairwise velocity dispersion u\2 measured at k = 
lh Mpc" 1 in the SDSS (circles, (3 = 0.45; squares, (3 varies with 
luminosity as in Tegmark et al. [2004]), compared with that in the 
2dFGRS (triangles) presented by JB04. The line without symbols 
shows the prediction of the HOD model of Yang et al. (2003). 



—21 < Mo.i r < — 20 sample) is constant around 500 kms -1 
at all scales. For faint galaxies, the PVD values increase as 
a function of k, reach a maximum value of 650 km s -1 at 
k ~ lh Mpc -1 and then decrease again. For bright galaxies, 
the behaviour is very different. The PVD values decrease 
as a function of k, reach a minimum value of 400 km s -1 
at k ~ 0.5h Mpc" 1 and then increase. It is interesting that 
these main features agree with the findings by JB04 for the 
2dFGRS, indicating that the PVD as a function of the scale 
k and the luminosity are robustly determined with these 
surveys. As shown by JB04, these features could not be re- 
produced by the halo model of Yang et al. (2003), because 
the model predicts that the PVD increases monotonically 
on small scales. Recently, Slosar, Seljak & Tasitsiomi (2005) 
made an attempt to interpret the observed luminosity de- 
pendence of the PVD in the context of a halo model, and 
showed that the luminosity dependence of the PVD can be 
predicted if some of the faint galaxies are the satellite galax- 
ies in high mass halos, which is consistent with what JB04 
expected as well as with galaxy formation models and HOD 
modeling (e.g., Guzik & Seljak 2002; Berlind et al. 2005; Ze- 
havi et al. 2005). However, a quantitative HOD model that 
can simultaneously match the observations, such as the lu- 
minosity functions, the correlation functions and the PVDs, 
still needs to be worked out. The observed features of the 
PVD are very sensitive to how the galaxies of different lu- 
minosity are distributed among dark matter halos and also 
inside these halos. The PVD is therefore a good constraint 
on both semi-analytical models of galaxy formation and the 
halo occupation models for the galaxy distribution. 

The luminosity dependence of the PVD is shown more 
clearly in Figure 6, where we have plotted o\2 at k = 
0.25, 0.5, 1 and Ah Mpc - for galaxies of different luminosi- 
ties. We again see that the observed PVD for M* galax- 



ies does not depend on fc-value. On large scales (k — 
0.25 h Mpc -1 ), <Ti2 rises as a function of increasing lumi- 
nosity. However, on small scales, the PVD exhibits a well- 
defined minimum at luminosities around M* — 1. The values 
of <Ti2 at this minimum are ~ 500 kms -1 at k = lh Mpc -1 
and ~ 400 kms -1 at k = 0.5 h Mpc -1 . We thus conclude 
that bright and faint galaxies have higher random motions 
than galaxies of intermediate luminosity. The maximum rel- 
ative velocities of faint galaxies (~ 600 kms -1 ) occur on 
scales ~ 1 h Mpc -1 . For the brightest galaxies, the maxi- 
mum values occur on the smallest scales and reach values 
close to 1000 kms -1 . 

It is interesting to compare the results presented in 
this paper with those of JB04. To do this, one must trans- 
form from the luminosities measured in the bj band in the 
2dFGRS to those measured in the 01 r band in the SDSS. 
According to the luminosity functions in the two surveys 
(Madgwick et al. 2002, Blanton et al. 2003), the average dif- 
ference in absolute magnitude in the two bands is ~ m .9. 
In Figure 7 we compare our SDSS results with those of the 
2dFGRS after taking into account this difference. It can be 
seen that the two measurements are quite consistent with 
each other. It is worth noting that the agreement is almost 
perfect for luminosities brighter than —19.5, but for fainter 
luminosities, <Ti2 is slightly smaller as measured from the 
SDSS compared to 2dFGRS. The latter may be due to the 
fact the volumes covered by the faint samples are not suffi- 
ciently large (Mo, Jing & Borner 1997). The results obtained 
with varying (5 are also plotted in this figure. We also plot 
the prediction for the 2dFGRS PVD of JB04 based on the 
halo model of Yang et al. (2003), which clearly does not 
match the observations. The CT12 measurements as a func- 
tion of luminosity are listed in Table 1. 
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Figure 8. The real space power spectrum P(k) (bottom panels; here the power spectra are given in the power-per-log-interval pre- 
sentation: A 2 (fc) = fe 3 P(k)/(2n 2 )) for galaxies in the luminosity interval of —21 < Mo.i r < —20 and with different properties (from 
left to right: g — r colour, D4000, concentration and log stellar surface mass density log 10 /Lt«), compared to the corresponding w p (r p ) 
measurements (top panels) from Paper I (see the third row in Fig. 10 of Paper I). In each panel, the black is for the full sample, while 
the red (blue) is for the subsample with larger (smaller) value of the corresponding physical parameter. 



3.3 The dependence on g — r, D40005 C and /u* 

To probe the dependence of the PVD on parameters re- 
lated to the recent star formation history of galaxies (g — r, 
D4000) and on parameters related to galaxy structure (C, 
jU*), we divide the galaxies into two subclasses: galaxies with 
larger values of the given physical quantity (hereafter de- 
noted as "red" galaxies) and galaxies with smaller values 
(hereafter denoted as "blue" galaxies). These subsamples 
were described in §2 of Paper I. 

Fig. 8 shows the real space power spectrum P(k) in the 
space of luminosity vs the physical quantities. The results are 
compared with the measurements of w p (r p ) presented in Pa- 
per I. As can be seen, the results of P(k) and w p (r p ) are qual- 
itatively very similar. Galaxies with redder colours, stronger 
4000A breaks, more concentrated structure and higher sur- 
face densities have higher clustering amplitude on all scales 
and at all luminosities, with the difference more pronounced 
on small scales and for faint galaxies. Furthermore, the de- 
pendence on parameters associated with recent star forma- 
tion is much stronger than the dependence on structural 
parameters. 

Fig. 9 shows the dependence of the PVD o\% on k 
for galaxies with different luminosities and physical prop- 
erties. The samples and the symbols are the same as in 



Fig. 8. In Fig. 10, we plot o\% measured at k =0.25, 0.5, 
1 and 4 h Mpc -1 as a function of luminosity for galax- 
ies with different physical properties. Plots corresponding 
to Figs. 9 and 10, as a function of stellar mass rather than 
luminosity are shown in Figs. 11 and 12. We find that red 
galaxy populations as defined by both g — r and D4000 
and early-type galaxy populations as defined by C and /i* 
have larger relative velocities on all scales and at all lu- 
minosities/masses than blue and late- type populations. We 
can also see that the dependence of 012 on these physical 
properties is stronger for faint/low mass galaxies and on 
small scales. On very large scales, galaxies of all luminosi- 
ties/masses have very similar pairwise peculiar velocities. 

In order to study the detailed dependence of the PVD 
on galaxy properties, we restrict our analysis to galaxies with 
stellar masses in the range 10 < log 10 M* < 11, and divide 
the sample into subsamples with finer bins in g — r, D4000, C 
and fi, (see Table 3 of Paper I). In Fig. 13 we show the pair- 
wise velocity dispersion CT12 on different scales, as a function 
of these quantities. The dependences of o\2 on quantities re- 
lated to past star formation history and on quantities related 
to galaxy structure are quite different. Fig. 13 shows that 
(T12 increases with increasing g — r or D4ooo- On the other 
hand, <7i2 as a function of concentration exhibits a maxi- 
mum at C ~ 2.6 for k = 0.25 h Mpc" 1 and at C ~ 2.8 for 
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Figure 9. The pairwise velocity dispersion <j\i for galaxies with various luminosities and physical properties. Panels from top to bottom 
correspond to galaxies in different luminosity intervals (as indicated in the first column), while panels from left to right correspond to 
galaxies divided by different physical quantities (as indicated above each column). In each panel, red (blue) is for the subsample with 
larger (smaller) value of the corresponding physical quantity, while black is for the full sample. 



larger k values, reaching ~ 700 km s -1 on small scales and 
~ 500 km s -1 on large scales. Furthermore, on very small 
scales, "star-forming" galaxies with D4000 < 1.5 have nearly 
constant peculiar velocities (~ 300 km s -1 ), while 012 in- 
creases sharply above -D4000 = 1.5, reaching ~ lOOOkms -1 
for galaxies with D4000 = 2. The PVD is an indicator of 
the depth of the local gravitational potential. Therefore, we 
are led to the conclusion that the reddest galaxies move in 
the strongest gravitational fields. A substantial fraction of 
these red galaxies must reside in clusters, while most galaxies 
with blue colours, recent star formation, and diffuse struc- 
ture populate the field. Interestingly, the behaviour of ai2 
as a function of concentration is more complicated. Galax- 
ies with intermediate value of C have the highest relative 
velocities. This may indicate that many red galaxies in clus- 
ters are really spirals whose ongoing star formation has been 
truncated. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we have made a detailed study of the power 
spectrum and the pairwise velocity dispersions of galaxies 
classified according to their luminosity, stellar mass, colour, 
4000A break, concentration index, and stellar surface mass 
density, using the SDSS DR2. Our results can be summa- 
rized as follows: 

(i) The real space power spectrum and the pairwise ve- 
locity dispersion as a function of luminosity in the SDSS 
are in good agreement with the result presented by Jing & 
Borner (2004) for the 2dFGRS, after taking into account the 
photometric bandpass difference between the two surveys. 

(ii) Galaxies with redder colours, stronger 4000A breaks, 
more concentrated structure and higher surface mass densi- 
ties have larger clustering power and larger relative velocities 
on all scales and at all luminosities/stellar masses. 

(iii) The dependences of the clustering power and PVD on 
the parameters related to recent star formation are stronger 
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Figure 10. The pairwise velocity dispersion 012 (fc) as a function of absolute magnitude, measured at k =0.25, 0.5, 1 and 4 h Mpc -1 
for different classes of galaxies. Panels from top to bottom correspond to measurements on different wavenumbers, while panels from left 
to right correspond to galaxies divided by different physical quantities, as incidated above each column. In each panel, red (blue) is for 
the subsample with larger (smaller) value of the corresponding physical quantity, while black is for the full sample. 



than those on the parameters related to galaxy structure, 
especially on small scales and for faint galaxies. 

(iv) The reddest galaxies and galaxies of intermediate 
concentrations move in the deepest gravitational field. A 
large fraction of these objects must reside in clusters, while 
most galaxies with blue colours, recent star formation, and 
diffuse structure populate the field. 

From the above results, it is not difficult to understand 
why the PVD on small scales (k = 1 h Mpc -1 ) exhibits a 
minimum at M* - 1, and why faint galaxies have larger rel- 
ative velocities than bright galaxies (JB04). As seen both 
in Fig.8 in this paper and in Fig. 10 of Paper I, faint red 
galaxies are more strongly clustered than faint blue galaxies, 



and have comparable or even greater clustering power than 
bright red galaxies. However, the majority of faint galax- 
ies have blue colours and cluster weakly. The fraction of 
red galaxies is only about 30 percent for luminosities fainter 
than —19 (Table 1 of Paper I). This means that the power 
spectrum of faint galaxies is dominated by the blue popu- 
lation, which leads to a monotonic increase of the cluster- 
ing strength with luminosity (Norberg et al. 2001; Zehavi 
et al. 2005; Tegmark et al. 2004). However, the faint red 
population in rich clusters, even though small in number, 
can dominate the PVD on small scales, because the PVD is 
more sensitive to the galaxy population in rich clusters than 
the power spectrum (see Mo, Jing & Borner [1997] for a dis- 
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Figure 11. The pairwise velocity dispersion 012 for galaxies with various stellar masses and physical properties. The symbols are the 
same as in Fig. 9, except that the green line in each panel is for the full sample volume-limited in stellar mass. 



cussion). This is also clearly seen in Figure 10. Our results 
therefore support the conjecture made by JB04 that a large 
fraction of the faint galaxy population must reside in rich 
clusters. We demonstrate that this fraction is predominantly 
red. 



Although most of the results presented in this paper can 
be understood qualitatively by considering known trends in 
galaxy properties as a function of environment, a quantita- 
tive understanding requires galaxy formation models (e.g. 
Kauffmann et al. 1993, 1997, 1999; Cole et al. 1994,2000; 
Somerville & Primack 1999; Kang et al. 2005; Croton et 
al. 2006) or halo occupation models (e.g. Jing, Mo, Borner 
1998; Seljak 2000; Berlind & Weinberg 2002; Cooray & 
Sheth 2002; Yang et al. 2003; Zheng 2004; Slosar, Seljak, 
Tasitsiomi 2006). Through such studies, the results in this 
paper will provide stringent constraints on these models, 
and will provide important clues where galaxies of different 
physical properties are located and how they have formed 
and evolved. 
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